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Abstract 

The s-wave interaction of DAci DT,c, D*Ac, D*T,c and DU*, D*T,*, is studied within a uni¬ 
tary coupled channels scheme with the extended local hidden gauge approach. In addition to the 
Weinberg-Tomozawa term, several additional diagrams via the pion-exchange are also taken into 
account as box potentials. Furthermore, in order to implement the full coupled channels calcula¬ 
tion, some of the box potentials which mix the vector-baryon and pseudoscalar-baryon sectors are 
extended to construct the effective transition potentials. As a result, we have observed six possible 
states in several angular momenta. Four of them correspond to two pairs of admixture states, two 
of DJlc ~ D*T,c with J = 1/2, and two of DT,* — D*T,* with J = 3/2. Moreover, we find a D*Jlc 
resonance which couples to the DAc channel and one spin degenerated bound state of D*Tf. with 
J = 1/2,5/2. 
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I. INTRODUCTION 


In this work, possible nucleon resonances with negative parity are studied, with the energy 
around 4200 ~ 4400 MeV, which are found as molecular states of one anti-charmed meson, 
D or D*, and a charmed baryon. Here the local hidden gauge approach with an extension 
to SU(4) is utilized, as done in Refs. to deal with the hadronic interaction involving 
heavy quarks. Within the framework of the local hidden gauge approach, the meson-baryon 


scattering proceeds via the vector-exchange 


(see also Ref. [8| for practical rules) and 


the interaction kernel projected over s-wave results in the form of the Weinberg-Tomozawa 
interaction 0, H Um- The heavy quark spin symmetry (HQSS) is known to play an 
important role in the heavy quark sector |l4Hl 



and has been utilized to investigate the 
1^2^. Yet, it was also found in Ref. [I| 


hadron spectroscopy in the high energy region [1| 


that the extension of the local hidden gauge approach automatically fulhls the rules of HQSS. 

Compared to the light quark sector, more importance should be placed on the vector- 
pseudoscalar admixture because of the small mass difference between them, namely mu* — 
rriD ^ 140 MeV in the charm sector while rriK* —rriK ~ 500 MeV. In the meson-baryon scat¬ 
tering, the mixing of vector-baryon {VB) and pseudoscalar-baryon {PB) sectors proceeds 
via the pion-exchange induced by the DD*7 t vertex, which is one of the ingredients of the 
local hidden gauge approach 0-0, 25l-27|. When dealing with pion exchange, the on-shell 
factorization of the potential, which is done for vector exchange in the chiral unitary ap¬ 
proach, is not possible because of the momentum dependence of the pion exchange. This is 
why the dynamics of pion exchange is included by means of a box diagrams, mediated by 



tion does not interfere with the s-wave driving force of the Weinberg-Tomozawa term, and 
hence it also goes via the pion-exchange in the box diagram. This box diagram involving the 
anomalous coupling contributes only to the VB sectors as an attractive potential. When a 
VB channel develops a bound state, the corresponding PB channel also generates a bound 
state and the generated state of the VB sector should be heavier than that of the PB sector 
due to the mass difference of D* and D. However, thanks to the extra attraction coming 
from the anomalous term, the VB sector binds more and hence the mixing effect of the VB 
and PB sectors now becomes more noteworthy. 

In this work, we go further with our approach to the VB and PB mixing by implementing 
the full coupled channels calculation. Indeed, the VB sector can couple to the PB sector via 
the pion exchange interaction, however this exchange carries momentum transfer and thus in 
the multiple scattering it is not straightforward to factorize the relevant terms concerning this 
VB-PB mixing. Then we implement the full coupled channels calculation by constructing 
effective VB ^ PB transition potentials from the corresponding box diagrams as done in 
Ref. 


Eventually, six nucleon resonances are generated within the framework of the full coupled 
channels calculation based on the extended local hidden gauge approach with the anomalous 
term with several quantum numbers. In addition to the masses and widths of the states, 
their building blocks or main components are evaluated by the wave function at the origin. 
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FIG. 1. Diagrammatic representation of the pseudoscalar-baryon interaction (a) and vector-baryon 
interaction (b). 


as discussed in Ref. 29|. One of the remarkable observations is the appearance of pairs 
of orthogonal states mixing the same PB and VB states. In the previous work on open 
charm baryons jl^, in addition to Ac(2592) which has experimentally been observed, its 
orthogonal state with an energy 2767 MeV was also generated. Similarly, in this work two 
pairs of orthogonal states are also deduced. Moreover, as HQSS is incorporated implicitly 
within the driving force, the results resemble pairs of heavy quark spin symmetry partners. 
The meson-baryon hidden charm sector was studied along the same line in Refs, j^, but 
without including pion exchange and hence not mixing the PB and VB states. We shall 
comment on the analogies and the differences found here with respect to that work. We 
shall also comment on results of other works that use SU(8) spin-flavour symmetry 
quark model 3 • 


23 


or 


II. FORMALISM 


A. The interaction via the vector-exchange 


We study the isospin J = 1/2 hidden charm states, the N* resonances with the ener^ 
around 4200 ~ 4400 MeV, following the local hidden gauge approach [sl-Q. In Refs. [I|, lil, 
the same states have been investigated within the local hidden gauge approach extended to 
SU(4). The validity of the use of SU(4) in this kind of studies has been widely discussed 


in Refs. [13|, l27|, l30|. Although SU(4) is broken at the end as a consequence of the multiple 


scattering and the different masses of the particles within a same multiplet, the elementary 
vertices which we use here are rather SU(4) symmetric and its use within the present context 
is acceptable. 

We will use the coupled channels formed with anti-charmed mesons D and D* and 
charmed baryons Ac, Sc, S*, Ac, namely DAc, DT,c, D*Ac, D*T,c, DT,*, D*T,*. In Ref. j^, 
the lighter states such as r]cN, ttA, rjN, rj N, 77S, KA for the PB sector, were included 
as a correction (a similar procedure was also considered for the VB sector), however, these 
states were found not to modify the energies and simply gave around 30 MeV width to the 
generated states. In the present work, we do not consider these states but will keep in mind 
that the states generated that can couple to these states will have a width of that order of 
magnitude at least. 

According to the local hidden gauge approach, the meson-baryon interaction proceeds via 
the exchange of vector mesons as depicted in Fig. [TJ Neglecting the three momenta carried 
by the exchanged vector, as is automatically done in these studies, we obtain the s-wave 
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potential of the PB scattering from channel i to channel j 31 


Vij — -Cij^j^{2y/s - 


Bi 


M. 


Bi 


MBi+Ei IMb,+E^ 
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( 1 ) 


with / = /tt = 93 MeV being the pion decay constant, Ei {Mb^ Ej) the mass, energy 
of baryon of i {j) channel. In the case of the VB scattering, the interaction kernels and 
the following scattering amplitudes have an extra factor e ■ e the scalar product of the 
polarization vectors of the incoming and outgoing vector mesons. The coefficients of Eq. ([T]) 
are given in Tables lAll IA21 1 A3 1 and IA4I in the Appendix. 

As we shall see in Sec. Ill Cl in addition to the vector-exchange interaction of Eq. ([T]), 
the pion exchange interaction is also taken into account which stems from the local hidden 
gauge scheme too. This contribution is responsible for mixing the PB and VB sectors and 
breaking the spin degeneracy as studied in Refs. (25-27, 32-34 . 


B. The construction of the scattering matrix 


With the interaction kernel Vij of Eq. ©. we evaluate the PB or VB two body scattering 
amplitudes. By solving the Bethe-Salpeter equation in coupled channels, we obtain the 
scattering amplitude T given in a matrix form 

T=[1-VG]-^V, (2) 


with G being the diagonal loop function which stands for the propagation of the intermediate 


meson-baryon states. As pointed out in Refs. [2|, Il3[, the use of the G function with 


the dimensional regularization scheme can generate unphysical states below the threshold. 
Hence, in this work, we utilize the G function regularized with the cut off scheme 


G{s) = 


^Qma. 


2M. 


B 


dP’q ojp -|- ojb 
(27r)3 2upUB P^"^ — {ojp + ubY + is' 


( 3 ) 


where is the CM energy, s = (P°)^, up = \/q‘^ + m?p, up = \/q‘^ + M|, and gmax is the 
cut-off of the three-momentum. 

Here, let us refer to a problem which occurs in the cut off method too. The G function 
with the cut off method also has deficiencies if the cut off is not reasonably larger than the 
on shell momentum of the intermediate state. This can only affect the DAc channel in our 
approach. But the interaction of the DAc channel is repulsive and does not lead to bound 
states in this channel. This channel has only a weak and indirect effect in our results which 
will be discussed below. 


C. Breaking the spin degeneracy in vector-baryon sectors 


In this subsection we break the spin degeneracy in the VB sector, namely that of the 


D*B states. Following the approach of Refs. [25|, l32|, |33|, |35 


we mix states of D*Ba and 
DBf) states in both spin channels as a box diagram such as D*Ba —)■ DBi, —)■ D*Ba and vice 
versa, in a first step. Later on we will revert to full coupled channels. The procedure to 
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FIG. 2. Diagrammatic representation of the in the intermediate state (left) and the DSc in 

the intermediate state (right). 


evaluate these diagrams follows closely what has been done before in the open charm case 


but in the baryon sector we have more baryons to take into account, while in Ref. 
only nucleons were involved. 

We consider the following PB ^ VB mixing processes 


27 


DSc -)• D*Ec -)• DEc, 

DEc D*Ac DT,c, 

DEI D*K DK, 

D*Ec -)• DAc D*Ec. 


D*Ec DEc -)■ D*Ec, 
D*Ac DEc -)■ D*Ac, 
D*El DEI D*EI, 


In the mixing we omit the transitions involving the vrScS* vertex. The reason is that the 
SU(4) Clebsch-Gordan coefficient for this vertex is l/\/l2 compared with l/\/2 for the vriVA. 
Given the moderate role played by vriVA transition found in Ref. [l^, this approximation is 
a fair simplifying option. The Glebsch-Gordan coefficient for TrA^S* is l/\/4, not so small 
compared to that of ttAA, but given the fact that the DA^ and D*Ac channels do not 
generate bound states by themselves and play only a minor, indirect role in the study, we 
also neglect terms involving this vertex from the beginning. In addition, note that the TrA^Ac 
vertex is zero due to isospin conservation. 

In Fig. [21 a pair of the DEc and D*Ec box diagrams, is shown. As discussed in Sec. IIIBl 
we include the DAc channel only as an intermediate state in the D*Ec box diagram. Then 
the channels with S* are completely separated from the others with Ec or A^. 

First we show the construction of the box potentials with the baryons, A^ or We 
detail the calculation for DEc D*Ec —)■ DEc and use the result as reference to obtain the 

result in the other cases. By using the D{D*) isospin doublet, {D^, D~) and Sc(F*) triplet, 
we write down the relevant hidden charm states with I = 1/2 in the charge 

basis 


5E„ / = 1/2,= +1/2) = + /1|C“E/-), (4) 

5A„ / = 1/2, /3 = +1/2) = |£)»A+), (5) 


and analogously for their excited states such as D*Ec, DE^, D*E*^ and D*Ac. In order to 
evaluate the DEc D*Ec transition in J = 1/2, we must consider the diagrams depicted 
in Fig. [3l and technically, as also done in Refs. |26|, |27(], we evaluate the transitions using 
only SU(3) symmetry by working in the {u,d,s) sector and using the correspondence of 
the and —K~) isospin doublets with the {D^,D~) and {D~^,—D^) ones. 
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FIG. 3. Diagrammatic representation of the transition DY^c X- D*Ylc in / = 1/2. 


respectively. We need the VPP Lagrangian 


Cvpp = -ig (|F,8„P]C'‘), 


( 6 ) 


where P, are the ordinary psendoscalar octet and vector nonet SU(3) matrices of the 
corresponding fields 


P 


V, 


72 V6 
n~ 

K- 

f A JiL 
72 ^ 72 

P~ 
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TT^ 

7r° m 
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p° I ^ 
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K*o 


K+ \ 

2r?8 
76 / 

i7*+ \ 

K*o 

^ / 


( 7 ) 

( 8 ) 


and (7 = my/^fT, with my 780 MeV. 

From the Lagrangian of Eq. ([S]) with Eqs. ([7]) and (|H]), the K* —)■ 77vr decay amplitnde that 
corresponds to the D* —)■ Du process is obtained. Yet, in order to get the conpling in the 
charm sector, it is reqnired to implement an extra factor to acconnt for the different meson 
held normalizations (l/\/2n7)and then fnlhl the HQSS rnles, as pointed ont in the previous 
study of the beauty sector |26[. In this case, for instance, the two transition amplitudes, 
D*~ —)■ D^n~ and K*^ —)■ K^ti~ ^ are related as follows 


^ ^ mp* 

tK*^K-K mx* 


( 9 ) 


The analogy of the strange and charm sectors used before for the ttD*D vertex is also 
used for the ttBB vertex, using the correspondence of the isospin multiplets (—S’*", S 7 ) 

and (—S°). Then we use the standard Yukawa coupling of pions to baryons in 
SU(3) [361, l37|, and we obtain the amplitude for the DTjc —)■ 7)*Sc transition 


iP 


= -V2g 


mp* 


= 2V2g 


rriK* 

mp. 


^ ^ 1 


1 


2F 


rriK*^ ''q‘^-ml2f, 


- ml 2 /, 

2E ^ ^ 


-a ■ q 


-a 


( 10 ) 


taking D = 0.75 and F = 0.51 [38j for the two couplings of the Yukawa vertex, where 
Pin, Pout are the incoming, outgoing meson momentum and q is the momentum transfer, 
and several relations. Pin = q + Pout and Pout -6 = 0 plus e° ~ 0, are used. Note that 
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FIG. 4. Diagram of the Kroll Ruderman term. 
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FIG. 5. The DY^c D*Y,c DY^c box diagrams. 


we work close to the threshold of the vector-baryon channel, which justihes taking = 0. 
In addition to the pion exchange of Fig. [31 we shonld also consider the Kroll-Rnderman 


contact term from the gange invariance constraint, as depicted in Fig . 
e^(g + Pin)^ g 2^^2 ^ ■ Q of Eq. (ITOll by —a ■ e, as discnssed in Refs. [25|, 
contact term 




/- mo* 2F ^ ^ 

= V2g - —a • e. 

tuk* 2/^ 


By replacing 
28| . we have the 


( 11 ) 


With the two DYc —)■ D*I]c transition amplitndes of Eqs. flTOll and ffTTl) . we constrnct the 
potential for the DYc —)■ £>*11^ —)■ box diagrams as shown in Fig. [5l and we obtain 


51/ = + 251/"^'-' + 51/ 


PC 


rCC 


( 12 ) 


where stands for the hrst diagram of Fig. [5l 2 5V for the two middle diagrams 

and for the last one. As shown in Refs. 111-0, we have the box potential of the 

DYc —)■ D*^c DYc process 


- iSV^P = 


d'^q (mo 
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2 uD*{q) - g° - ujD*{q) + ^s.(g) - ^s.(g) + ie 


, ( 15 ) 


where Pj^ and are the energies of the incoming meson D and baryon E,, respectively, 
and uJoY^ = ^/q^ + ml,, etc., Ps,(<h) = sjq"^ + m\^. 

It is easy to handle the operators a and e in the two pion exchange term of Eq. for 
both cases of the PB —)■ EP —)■ PB and EP —?■ PB VB boxes. In the hrst case we can 
nse the identity {q ■ = q^, and then the symmetric property in the integral 


QiQj 

and the snmmation of the polarization of the intermediate vector 


(16) 


= iii- 


(17) 


However, in the HP —?• PB —)■ HP box, the polarization vectors are external and there is 
no snm over polarizations, one is e and the other e'. We still can nse Eq. flT6|) and we obtain 
in both cases 


{{a-q){e-q)V^ 


1 q^e- e' 


PB^VB^ PB 
VB ^ PB ^VB • 


(18) 


From this resnlt, the inverse box diagram, D*'Ec —)■ PSc —)■ P*Ec, is fonnd to have an extra 
factor 1?• e'. Furthermore, for the contributions which include the Kroll-Ruderman term. 


namely Eqs. (na and fllSp . we use again Eq. fll6p plus the property [23 


(FB|c?-e1l7S> = VS6j,i/2. 


(19) 


and we hnd in both cases that there is only J = 1/2 contribution and 

{3 ■ e){a ■ ■ if) =q'^h,\l2- 


( 20 ) 


For other box diagrams involving t-)- baryon conversion, one can obtain the potential 
in the same manner as done above. The hnal result is that we can use Eqs. (na-ci, 
considering the appropriate masses or energies and making the following replacement 
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( 21 ) 


Next, we evaluate the box diagrams of the sectors involving a baryon S*. The two 
diagrams, PSc D*Jlc —>■ PSc and PS* —)■ P*S* —)■ PS*, are different in their vrPP 
vertices. The 7rS*S* vertex is given by 


- itnT)*E* 


/s* 




Ss* • q Ts* 


( 22 ) 


where Sy,* and Ts* represent a spin S' = 3/2 and isospin 1 = 1 operator respectively. Using 
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the ordinary quark model, one finds 


/sg _ 4 fnNN _ 4: D + F 

5 5 2/ 

By counting the weight of each contribution in the same manner as for the case of Sc, we 
have the -DS* —)■ -D^S* transition amplitude 


it 




2 V2g 


TflD* 

-g • 

rriK* 



— ml 



(24) 


A further simplification is done in this case and we evaluate the box for an average of the 
diagonal spin transition. We have 

^ ■ ^m) = '^{m'\S^*\m)i{m\S^*\m')j. (25) 

/ / 

Since there is no privileged direction after we sum over m,m' the matrix element, we have 

^ ^ (m'|S's«|m)i(m|S'E.|m')j = A5ij, (26) 

/ 

m,m 


and by taking the trace, we obtain the spin factor A = 5/4. This allows us to write the 
following relationship 

■ gj ■ g) -t -g^. (27) 

Recalling the identity, {a ■ q)'^ = q'^ used in the case of spin 1/2 baryons in the intermediate 
states, in order to construct the box potential of the Z)S* —)■ Z)*S* —)■ /)S*, we implement 
the following substitution in Eqs. ffT^ - flTSl) 


2F 

V 




(28) 


D. Contribution from the anomalous term 

In addition to the VB and PB mixing term, another correction contributes to the D*B 
channels which does not interfere with the s-wave driving force at tree level. This contribu¬ 
tion stems from the anomalous D*D*tt coupling and is taken into account in a box diagram 
too. By replacing the D mesons in the PB-VB box by D* ^ we list down the possible cases 
for this anomalous box contribution 

T)*Sc ^ i)*Sc ^ i)*Sc, 

Zl*Sc ^ D*Ec, 

D*Ac D*Ac, 

-)■ -)■ 

One of them, the 7)* Sc —)■ 7)*Sc —)■ 7)*Ec diagram, is depicted in Fig. [6l 
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FIG. 6 . Box diagram with anomalous D*D*tt vertex. 


The anomalous D*D*tt vertex can be obtained from the Lagrangian0 


Cvvp = 

v2 


(29) 


with G' = 3My/167r^/^(~ 14 GeV j39|, l40[ and the V and P matrices in SU(4) are given 
in Ref. M|- Thus, for the D*E* — D*E* transition amplitude we obtain 


1 G' 2F 1 

“'^^D*Sc^-S*Sc ^ - 2 \/ 2 ^m 25 *e ei{p} qj _ ^2 


a ■ q, 


(30) 


with a four momentum of the incoming Z)*, ~ ijriD*) 0), and = e^“^. 

In the case of the anomalous term, there is no contact term. Here the box potential is 
obtained in analogy with the case of the previous box diagram, DY^c —)■ —)■ DYc- By 

using Eqs. flTB]) and flTTj) . the operators of the intermediate states contract and lead to 

" = pf. r' - («-■ £-)(«. £-') ^ (31) 

This, together with = q‘^ leads to a global factor ^q^e-e' instead of the factor ^q‘^e-e' 

that we had in Eq. (flSj) for the case of intermediate vectors. 


III. FULL COUPLED CHANNELS CALCULATION 

With an extension of the mixing effect between PB and VB channels caused by the pion 
exchange, we implement the full coupled channels calculation. Recalling the box potentials 
constructed in Sec. Ill Cl we finally study the following four sectors separately. 

J = 1/2 : DY^, D*Y^, D*Ac 

J = 3/2 : D*Y^, D*A^ 

J = 1/2 : D*Yl 
J = 3/2 : DYl, D*Yl 

The driving force, the Weinberg-Tomozawa term, is already given in Sec. Ill Al and it is found 
that these interactions have no off-diagonal interaction (see Tables lATl and lA^ . Furthermore 

^ An anomalous process, like the VVP interaction, is one that does not conserve “natural” parity. The 
“natural” parity of a particle is defined as follows: it is +1 if the particle transforms as a true tensor of 
that rank, and —1 if it transforms as a pseudotensor, e.g. tt, 7 , p and ai have “natural” parity — 1 , -|- 1 , +1 
and — 1 , respectively. 
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the other interactions, the box diagrams from the local hidden gauge and the anomalous 
term, are added. Since a large mixture oiVB and PB states would be expected from the re¬ 
sults of Ref. 27| , it is preferred to implement the full coupled channels calculation, especially 
between VB and PB channels which generate physical states or have attractive interaction. 
In this section, hrst we show the derivation of the effective transition potentials which allow 
us to implement the full coupled channels calculation, and subsequently the explicit form of 
the additional interactions to Weinberg-Tomozawa one are given. In Tables IA61 IXTI and 
lASl of Appendix, the interactions for each sector are summarized. 


A. Effective transition potential 


We follow the scheme of the effective transition potential constructed in Ref. j27[. In 
this work, the generated states are looked for in the s-wave scattering amplitude, where the 
Weinberg-Tomozawa interaction appears. However, the intermediate states of the box po¬ 
tentials are not necessarily in the s-wave and then the box potentials should be decomposed 
into two parts to implement full coupled channels scheme. The part with intermediate states 
in the s-wave is utilized to construct the s-wave effective transition potential and the rest 
part is added to the Weinberg-Tomozawa interactions as a correction. On the other hand, 
since the anomalous D*D*tt couping does not interfere with the Weinberg-Tomozawa term, 
these terms contribute only via the box diagram. 

With the s-wave component of the PBa VB^ ^ PBa box potential, 6 VpB^^vB^^PBa{^~ 
wave), such as DUc Z)*Ac(s-wave) —)■ we dehne an effective transition potential 

between PB„ and VBh channels as 


- ^^PBa 


> VS,(s-wave). 


(32) 


with GvBi being the G function of the intermediate VBi, channel. In the same manner, from 
the counterpart diagram, VBf, ^ PBa —t VBb, we can also dehne the other PBa ^ VB), 
transition potential with a corresponding box potential 6 V 


vliGpB^ = 6r 


VBh^PB, 


,-,VB,(s-wave). 


(33) 


To ensure the symmetric property of the potential or the amplitude matrices, namely Vij = 
Vji, we dehne the ehective transition potential by taking the average of the two potentials 


Vff - 2(Vff + Kff)- 


(34) 


However, it turns out that the potentials are very similar and in fact they are identical 


in the limit oi mp) = mp)* as demonstrated in Refs. [2^ l27|, hence, the concept of Ws is 


well dehned. From the dehnitions, We and are found to be double-valued functions. 
Therefore, the relative signs of the two functions should be taken equal not to be canceled 
out. We choose I4fr with negative real part as it corresponds to a virtual pion exchange. 
The energy of the generated states does not depend upon this prescription. 

Among the relevant box diagrams, it is obvious that the term involving the contact term, 
Eqs. flTT)) and flT^ . contribute only to the s-wave due the product a ■ e. Therefore, the 
intermediate states of these terms are always in the s-wave. On the other hand, the case of 
the two pion exchange term is different. When the extraction of the s-wave component is 
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necessary, we utilize the following separation for each of the pion exchange 

(? ■ 9 ) (e ■ 9 ) -» ei9i9j9j = I ’ 


( 36 ) 


where the hrst term of the last expression of the equation corresponds to the s-wave com¬ 
ponent and the second to d-wave, and this separation leads to 


{(d ■ g) (e ■ g)}^ (s-wave) ^ <j 11 hL ^ 

3 3 y ^ ^ 


; PB ^VB ^ PB 
' : VB ^PB ^ VB 


(36) 


Subtracting the s-wave contribution from the total contribution of the box we obtain the 
d-wave contribution. 

In the following sections, the explicit forms of the potentials are shown where the ex¬ 
pressions of I 2 , 13 , FAC, Ii, I 2 , and I 3 can be seen in Eqs. (34)-(37) and (41)-(43) of 

Ref. 0. 


1. DTic, D*T,c, D* Ac with J = 1/2 

In this sector, the VB intermediate states in the PB —)■ RR —)■ PB box potentials can 
have s-wave and d-wave contributions. By using the decomposition, Eq. fl36l) . and recalling 
that the contact-contact and contact-pion terms contribute to only the s-wave component, 
we can decompose them into two parts 


dV/(s-wave) = 6 V (DSc —)• D*'Ec —>■ DT^c; J = 1/2, s-wave) 

( 1 

-^ + 2/2 + 13 
3 

dV//(s-wave) = 5V (ZdEc —)■ D*Ac —)■ DTjc, J = 1/2, s-wave) 

= REL2 X FAC x ( + 2/2 + h 

3 omi 


(37) 


(38) 


5V {DFc —>■ D*'Ec DFc, J = 1/2, d-wave) = RE LI x FAC x - 

6 V (DEc D*Ac DEc, J = 1/2, d-wave) = REL2 x FAC x „ „ , 

^ ^ 3 omi 


2 dh 

3 ’ 
2 dh 


(39) 

(40) 


where we use the factor FAC together with the relative factors as discussed in Sec. Ill Cl 

'F + D 


2 I '^D* 


FAC = -g^ 

2 \mK* 


2 / 


RELl = - 


2F 


REL2 = 


9 \D + F 


2D 


(41) 

(42) 

(43) 


9 \D + F^ 

The PB intermediate states in the VB ^ PB VB box potentials are necessarily in 
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the s-wave and thus with Eqs. fflSj) . fH2|) and fH3|) we have 

SVj{s-wave) = dV {D*E^ DS, ^ J = 1/2) 

= RELl X FAC X + 2/; + Z^) , (44) 

6 Vjj{s-wave) = 6V {D*A^ DS, ^ D*A^-, J = 1/2) 

= REL2 X FAC x + 2/' + . (45) 


From the s-wave box potentials given by Eqs. fl371) and fj441) and by Eqs. 
we construct the following two effective transition potentials 


^ DT,cAAD*Y.a 


^ DT,cAAD*Ka 


Y/(51//(s-wave)/Gp*s, + \/^i (s-wave)/Gf) 




\/ SVii{ s-wave)/Gp*yVc + V (s-wave) / G5 




and (05]), 

(46) 

(47) 


They are placed in the corresponding off-diagonal channels. 


Furthermore, as discussed before, the DA^ contribution is just added to the D*Tic channel 
as a box correction. 


5V (T)*Se ^ DA^ J = 1/2) 


REL2 X FAC x 


\dml 


A 2Jn -1- Jc 


(48) 


We also show the box potentials which stem from the anomalous term 

dl' 

(514n ^ ^ D*J:c) = RELl x AFAC x (49) 

dl' 

514n ^ ^ T)*Se) = REL2 x AFAG x (50) 

dl' 

(5Kn (A>*Ac ^ ^ D*AA = REL2 x AFAG x (51) 

with a factor which comes from Eqs. fl30p and fl3ip 

AFAC = ^G'^ ml.. (52) 

These box potentials are added to the 7)*Ec and D*Ac diagonal Weinberg-Tomozawa inter¬ 
action equally in the J = 1/2 and J = 3/2 sectors. 
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2. I)*Ac with J = 3/2 

Since the PB intermediate states with J = 3/2 are automatically in the d-wave, we just 
add the following box diagram potentials to the diagonal interactions 

dl' 

SV (D*Sc ^ DE, D*Sc; J = 3/2) = RELl x FAC x (53) 

dl' 

5V (D*A, DEc ^ D*Ad J = 3/2) = REL2 x FAC x 7 ^., (54) 

dl' 

6 V (D*Sc ^ E>Ac ^ E>*Ec; J = 3/2) = i?EL 2 x FAC x 7 ^. (55) 

om^ 

Remember that the box potentials from the anomalous term, Eqs. fl49l) - fl5Tl) . are also added 
to this sector. Note that in this sector there is no transition between the F*Ec and D*Ac 
channels and hence, in this sector we do not have coupled channels but two single channels. 


3. with J = \I2, hjl 


In this sector, we have only a channel. Considering that the FE* channel with 

J = 1/2, 5/2 can not be in the s-wave, we hnd that the box potentials should include 
only the pion-pion term with Eq. fllSp . Taking into account Eq. fl2SD . we have the Z)*E/ —)• 
FE/ —)■ F*E* box diagram 


with 


5V (F*E* ^ FE* ^ F*E*; J = 1/2) 


REL3 X FAC x 



REL3 = 




D + F 


16 

45‘ 


(56) 


(57) 


This is just added to the D*YA^ diagonal Weinberg-Tomozawa interaction with the extra 
contribution from the anomalous term 


(f*e: ^ f*e: ^ f*e:) 


REL3 X AFAC x 



(58) 


4 . FS* and F*S* with J = 3/2 

We here consider the FE* and F*E/ coupled channels. The FE* intermediate state with 
J = 3/2 has an admixture of the s-wave and d-wave components and thus each potential is 
separated into s- and d-wave contributions as done in the PB-VB sector with J = 1/2 

6 Vijj{s-wave) = SV (F*E* ^ FE* D*El; J = 3/2, s-wave) 

( 1 o 

3 dm^ 


+ 2/0 -l- F 


(59) 
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and 


2 dr 

6 V -r DTZ -r J = 3/2, d-wave) = REL?> x FAC x -7r4- (60) 

On the other hand, althongh the Z)*S* intermediate state of the -DS* —>■ D*YA^ —)■ -DS* box 
diagram also has the s-wave component, the many spin contribntions of complicate 

the analysis. Therefore, we adopt only the dehnition of Eq. (I3^ instead of Eq. fl3T)) to 
constrnct the effective transition potential 

“ y^5Ejjj(s-wave)/Gps*5 (61) 

and by nsing this effective potential, we extract the s-wave contribntion from the total 
DS* ^ DEI box 

5V [DE* —>■ D*E* —>■ DE*; J = 3/2, except s-wave) 

= REL3 X FAC x + 21, + /a) - (62) 

We also add the potential from the anomalons term, Eq. fl58|) . to the D*E* channel. 


B. Resonance measured on the real axis 

In the fnll conpled channels scheme, we wonld see the generation of several physical states. 
Under the conpled channels nnitary scheme, the corresponding poles are looked for in the 
second Riemann sheet in the complex plane. In the case of bonnd states, poles appear on 
the real axis below the threshold in the hrst Riemann sheet and thns we look for the bonnd 
states poles as nsnal. From the residne of the scattering amplitnde at the pole, Mr, we 
obtain the conpling constant Qi, gj 

giQj = lim {y/s - MR)Tij{3/s). (63) 

With the conpling constant and the G fnnction at the bonnd state pole, we evalnate the 
wave fnnction at the origin of channel i as giGi^Mji) jj^ . 

Resonance poles, on the other hand, emerge in the complex energy plane. As we do not 
know how to extend the box potentials or effective potentials to the complex energy plane, 
we analyze the resonance states in another way. Instead of looking for a pole position of the 
scattering amplitnde, we search a peak position Mr on the real axis, where the sqnare of 
the magnitnde of the amplitnde |T(i/i)p has a maximnm. Since this position is not exactly 
the same between different scattering channels, we determine it in the dominant channel. 
Together with the peak position, we measnre the width of the state, Tr, on the real axis as 
being a distance between two points where |T(y^)p has a half maximnm valne. 

Next, we determine the conpling constants of resonance states with the measured mass 
and width. On the real axis and in the vicinity of the peak position, we assume the Breit- 
Wigner amplitude 
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At the peak position, the coupling constants of the dominant channel d is given by the 
imaginary part of the amplitude 


ImTjj(M^) 


9i9j 

rV2 


) 9d 



ImTdd(A/if) 


(65) 


where gd is taken as a positive number. In order not to lose the relative signs, the coupling 
constants of the other channels are determined from the ratio of the imaginary part of the 
amplitude in different channels 


_ lmTid{Mii) 

~ ImTddiMnf’^ 


( 66 ) 


The wave function at the origin, or the probability, can be obtained from the coupling 
constant and the G function at the pole position. In this case, since we determine the 
coupling at the peak position, for the sake of consistency we evaluate the wave function at 
the origin of channel i as giGi{Mii). 


IV. NUMERICAL RESULTS OF THE FULL COUPLED CHANNELS CALCULA¬ 
TION 

In this section, we show the results of the full coupled channels approach. In order to 
obtain the numerical results, we take three cut off parameters, in the pseudoscalar- 
baryon loop, in the vector-baryon loop and q^^ in the box diagrams. Since no nucleon 
resonance state in the high energy region of our interest has been reported, we tentatively 
adopt the cut off parameters used in the previous work of open charm baryons [1^ where the 
parameters were chosen to reproduce two resonances, Ac(2592) and Ac(2625). In Table [U 
three sets of three cut off parameters are listed. Furthermore, as well as in the previous 
work, we put the Yukawa form factor, F{q) = h?/{h? + on the ttEB' vertices in all the 
box diagrams, with cut off A = 1000 MeV, a standard value. In the following subsections, 
we show the results of the and D*T,c states and and states separately. 



set I 

set II 

set III 

Qmax 

600 

800 

1000 

'imax 

771 

737 

715 

9max 

527 

500 

483 


TABLE I. Cut off parameters, gmax(9max) the vector-baryon (pseudoscalar-baryon) G function, 
9max the box potentials. 


A. DSc and states 

In Table [TTl the peak positions and the widths of the resonant states with J = 1/2, 
generated by the coupled channels of DUc, and D*Ac with DA^ box diagrams, are 
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listed. As discussed in Sec. IIII Bl in this work, we extract the nature of the resonant states 
not by looking for the corresponding poles in the complex energy plane but by measuring 
the scattering amplitude on the real axis. We determine the peak position and width of 
generated states from the D*T,f. —)■ D*T,c scattering amplitude. 

We obtain two different narrow resonances. With the ambiguities from the choice of the 
cut off parameters, the lower energy state appears around 4228 MeV ± 15 MeV with about 
20 MeV width, while the higher energy state is very close to the DAf. threshold, 4295 MeV 
with 10 MeV width. The coupling constants Qi and the wave functions at the origin giGi 
for the middle set of the cut off parameters, obtained from Eqs. fl65l) and fl66l) . are listed in 
Table uni From these wave functions, it is found that two states appear as two orthogonal 
states roughly as [DJ^c ± D*T,c). This situation is similar to the one observed in Ref. 

where the -^{DN — D*N) corresponded to the Ac(2595), while there was a prediction for 
an orthogonal state about -^{DN + D*N) which was much less bound E|. This feature was 


also observed in Ref. [2l| where the higher energy state appears with zero binding 

It is also remarkable that for both orthogonal states, the D*Tjc component has a bit 
stronger coupling than the D'Ec component. This could be one of the consequence of the 
heavy quark limit. Although, even in the charm sector, the mass difference between D and 
D*, rriD* — mo ~ 140 MeV is not so small yet, the vector-baryon system binds more than the 
pseudoscalar-baryon due to the anomalous term contributions and this mechanism makes it 
possible that the vector-baryon and pseudoscalar-baryon systems couple easier. 



set I 

set II 

set III 

peak 1 

4241.7 

4227.6 

4218.6 

width 1 

19.5 

21.1 

21.5 

peak 2 

4296.8 

4295.1 

4294.5 

width 2 

13.1 

10.6 

9.6 


TABLE II. Peak positions and their widths in the full coupled channels of 
Z)Sc[4321], i)*Sc[4462], Z)*Ac[4295] with J = 1/2 measured in the D*Tjc —)• D*Tic channel 
as a function of (The number in brackets after the channel indicates the threshold mass of 

the channel. Units: MeV) 


(4227.6,21.1) 

DE, 


D*A, 

9i 

4.40 

5.39 

0.39 

Qi Gi 

-15.66 

-24.17 

-3.31 

(4295.1,10.6) 



D*A, 

9i 

-1.27 

2.28 

-0.11 

QiGi 

8.46 

-12.60 2.09 + i0.04 


TABLE III. The coupling constants to various channels for the poles in the J = 1/2 sector of D*Tic 
and coupled channels, taking the cutoff Set II of Table [H 


2 In Ref. 


27| the sign of VeS was taken positive and here we take it negative. The energies do not depend 


upon this prescription but the ± D*N) combination would correspond now to -^{DN =F D*N). 
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Next we see the D*'Ec and D*Ac single channels with J = 3/2. As expected, the D*Ac 
channel generates no state and we observe the D*'Ec broad resonance whose mass is 4217 
MeV and width is aronnd 100 MeV as listed in Table HVl This resonance is generated 
below the DTj* threshold bnt can decay into DAc via the pion exchange, as well as the two 
resonance states with J = 1/2. Compared to J = 1/2 resonances, this J = 3/2 resonance is 
mnch broader becanse of the difference of the D*'Ec —)■ DAf. D*'Ec box diagram as seen 
in Eqs. fl4S]) and dSB]) . 

Apart from masses or widths of the generated states, it is remarkable the similarity 
between the results of the hidden and open charm baryon. The vector-baryon system whose 
interaction is attractive, such as in D*Ilc or D*N, generate deeply bound states with J = 1/2 
and J = 3/2 and the pion exchange removes their spin degeneracies. On the other hand, 
the corresponding pseudoscalar-baryon system such as or DN generates shallow-bound 
states with J = 1/2. As a consequence, in the J = 1/2 sector, a pair of orthogonal states 
emerges as a mixture oiVB and PB bound states, while in the J = 3/2 sector one deeply 
(quasi-)bound state is generated. 



set I 

set II 

set III 

peak 

width 

4250.5 

140.8 

4217.7 

103.2 

4205.8 

82.0 


TABLE IV. Peak positions and their widths in the Sc [4462] single channel with J = 3/2 sector 

as a fnnction of qm’ax ■ (Units: MeV) 

We come back here to the results of Refs, jl, and admit an extra width of about 30 
MeV for the two lower -DSc) states from the coupling to the light PB or VB 

sectors. 

B. DT,* and D*Y,* states 

In this subsection we show the results of the generated states involving DTj* and D*I1* 
channels. In these sectors, bound states are developed and then we look for the corresponding 
poles as discussed in Sec. IIIIB1 In Table IVj the pole positions of states generated by the 
single channel with J = 1/2, 5/2 are listed. We obtain a bound state with an 

energy 4344 with 10 MeV uncertainty from the cutoff setup, which appears around 180 MeV 
below the threshold. 



set I 

set II 

set III 

Pole 

4354.5 

4344.1 

4337.5 


TABLE V. Poles in the L)*Sc[4527] single channel with J = 1/2, 5/2 as a function of cut off 
(Units: MeV) 

Next we see the results of coupled channels of and with J = 3/2. The 

pole positions of the generated states are listed in Table IVTl In this sector, we obtain two 
bound states which are not far from the bound states with J = 1/2. The two states 
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are separated by about 50 MeV and the higher state appears a few MeV below the DS* 
threshold. We show the coupling constants and the wave functions at the origin with the 
cut off set II in Table IVIIl From these wave functions, it is found that two bound states are 
the mixture of and D*T,* corresponding roughly to the ±Z)*S*) combinations. 



set I 

set II 

set III 

Pole 1 

Pole 2 

4330.6 

4384.1 

4324.9 

4377.8 

4319.9 

4374.4 


TABLE VI. Poles in coupled channels of L)E*[4385], i4*S*[4527] with J = 3/2 as a function of 
(Units: MeV) 


4324.85 

DK 


9i 

3.61 

4.89 

3,G',' 

-16.54 

-24.22 

4378.84 

DK 


9i 

-1.25 

3.00 

9.Gl^ 

12.30 

-17.57 


TABLE VII. The coupling constants to various channels for the poles in the J = 3/2 sector of 
DT,*, taking the cutoff Set II of Table HI 


V. SUMMARY AND CONCLUSIONS 


In this work, the dynamics of an anti-charmed meson and a charmed baryon is stud¬ 
ied in order to investigate the hidden charm resonances with isospin I = 1/2, N*, in the 
energy region around 4200 ~ 4400 MeV. Since all the hadronic particles of our interest 
contain the charm quarkjthe extended local hidden gauge approach with SU(4) symmetry, 
as studied in Refs. [H, IH, [^ , is employed to extract the information of the relevant vertices. 
Consequently, the leading interaction proceeds via the vector exchange and results in the 
Weinberg-Tomozawa interaction. Furthermore, we also take into account two types of addi¬ 
tional interactions, as a box diagram correction. One of them with the DD*tt vertex from 
the local hidden gauge scheme, proceeds via the pion exchange and requires the correspond¬ 
ing contact term. The other one comes from the anomalous D*D*n interaction. The former 
additional interaction allows the PB and VB sectors to couple to each other. As pointed out 
in the previous works, we also place high importance on this mixing effect. Hence we have 
solved the Bethe-Salpeter equation under the full coupled channels scheme by constructing 


the effective transition potential following the discussion of Ref. [27| . The particular values 


of the ttBB vertices, allowed us to classify the states into four sectors of the hidden charm 
systems, DT,c, D*T,*, D*Ac with J = 1/2, D*T,c, D*Al with J = 3/2, D*YA^ with J = 1/2 
and D*T.* with J = 3/2. 

By looking for the bound state poles and measuring the resonant peaks on the real axis, 
as discussed in Sec. IIIIBl we have seen that six states appear with several angular momenta. 
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Furthermore, from the wave function at the origin, the dominant components of the states 
have been evaluated. The properties of the states, masses, widths, dominant components 
and main decay channel, are summarized in Table IVIIIl 


main channel 

J 

(E, r) [MeV] 

main decay channels 


1/2 

4228, 21(51) 

DAc 


1/2 

4295, 11(41) 

DA, 

D*T.^ 

3/2 

4218, 103 

DA, 

D*T.l 

1/2, 5/2 

4344, 0 

— 


3/2 

4325, 0 

- 


3/2 

4378, 0 

- 


TABLE VIII. Energies and widths of the obtained states with the dominant component and main 
decay channels of each state. All the states are nucleon resonances with negative parity and have 
an estimated uncertainty of ±20 MeV. The numbers in brackets for the first two states correspond 
to the estimated width, adding the 30 MeV width obtained in Refs. BB from coupling to the 
light PB or VB sectors. 


The three generated states with Sc are, a pair of orthogonal states with J = 1/2, 
^(Z)*Sc ± DT,c), and, a resonance dominated by the D*T,c component with J = 3/2. 
The two orthogonal states with energies around 4228 and 4295 MeV respectively, differ 
from each other by about 65 MeV and have a small width due to the weak coupling to the 
Z)Ac. On the other hand, the l)*Sc state with J = 3/2 emerges around 4218 MeV with 
a large width 103 MeV. It is instructive to see the similarity in the generated states with 
S*. Around 100 MeV above the region of the three state of Sc, a pair of the orthogonal 
states with J = 3/2, ;^(Z)*S* ± l/S*) and one spin degenerate state with J = 1/2, 5/2 of 

DS* develop. Considering the mass difference between Sc and S*, Ms* — Ms^ ~ 65 MeV, 
this result might indicate the realization of the heavy quark spin symmetry regarding the 
charmed baryon pair. Sc and S*. 

In the end, let us discuss the result comparing them to those of Refs. 0,1,1. Apart from 
the use of the different cutoff parameters and the neglect of the light hadron sectors that 
simply give the width to the generated states, we employ the same interactions for the PB 
and VB sectors respectively. In Ref. 1, the PB sectors with the same quantum number 
develop physical states around 4250 MeV, not so far away from our results, however, the 
VB states appear much higher than ours, about 200 MeV higher. This is also the case in 
the work done in Ref. 4^ done along the same lines. 

In all these works 11,1 , the box diagrams were not considered and the PB and VB 
sectors were studied independently. We have seen in the present work that due to the pion 
exchange that induces PB and VB mixing, the VB states get much bound. In addition, the 
pion exchange in the anomalous sector, which is only operative for the VB sector, provides 
a further contribution to the VB binding. 

Another work that deserves attention in this discussion is the one of Ref. 
unitarized coupled channel approach, also mixing PB and VB states is done. The authors 


23 . There a 
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use the Weinberg-Tomozawa interaction as leading contribution and an extended SU(8) 
spin-flavour symmetry constraint, with modiflcations implemented to respect HQSS. They 
obtain also bound states but the energies are almost 200 MeV smaller than what we And in 
the present work. 


There are other calculations based on the quark model 


The authors also And bound 
states for different spins, and the energies range between 3900 ~ 4500 MeV. One should 
note that in this latter work, the predictions differ by about, 200 ~ 300 MeV depending on 
the interaction used, a color magnetic interaction in one case and a chiral interaction in the 
other. A detailed discussion of these different model is done in Ref. 
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The comparison of the results obtained here with those of other works is instructive and 
indicate the large uncertainties in the predictions in this sector. In this respect, one should 
mention an earlier work where states are found with a binding of 1000 MeV which. 


as discussed in Ref. j^, is somewhat extreme. These large theoretical differences contrast 
with the situation in other sectors, like the one of open charm, where the constraints of 
reproducing existing experimental data makes the predictions of different approaches more 
similar. Yet, one thing is shared by all these approaches and this is the existence of many 
bound states of the N* type. 


The former discussion should be sufficient incentive to search for these states, which 
would bring a clear evidence of exotic baryonic states, with a quark structure significantly 
different from the standard one of the baryons made up of three quarks. 
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APPENDIX: SUMMARY OF THE INTERACTION 


In this appendix we write the Cij coefficients of Eq. ([T]) for the different channels. 


C^J 

DJlc 

US, 

1 


TABLE Al. Cij coefficients for DT^c with 1 = 1/2 and = 1/2 . 
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L)*Se 

D*A, 

D*Sc 

I 

0 

D*A, 


-I 


TABLE A2. Cij coefficients for D*Tjc and D*Ac with 1=1/2 and = 1/2 ,3/2 . 



DK 

DS* 

1 


TABLE A3. Cij coefficients for US* with 1=1/2 and = 3/2 . 



D*S* 

D*T,l 

1 


TABLE A4. Cij coefficients for D*T,* with 1 = 1/2 and = 1/2 ,3/2 ,5/2 . 



D^c 

D*Se 

D*A, 

DJlc 

Eqs. { O + (I39D + (gOD } 

Eq. dll]) 

Eq. ii) 

D*^c 


Eqs. { (P) + is]) + dM]) + dsni) } 

0 

D*A, 



Eqs. {dU) + dSH) } 


TABLE A5. Interactions in the coupled channels DTjc, D*Tic, D *with = 1/2 . 



D*T.c 

D*Ac 

d*j:c 

D*K 

Eqs. {(P + dSl + dSSIl + dMD + dSQl)} 

0 

Eqs. { dm + dSD + dSID } 


TABLE A6. Interactions in the two single channels Il*Sc,Il*Ac with = 3/2 . 



D*S* 

D*jj* 

Eqs. { dU + dSSj) + dSSI) } 


TABLE A7. Interactions in the single channel D*T,* with = 1/2 ,5/2 . 
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DK 


DK 

Eqs. { (P) + ([62]) } 

Eq. ([H]) 

D*T,* 


Eqs. { O + iO]) + (IMI) } 


TABLE A8. Interactions in the coupled channels with = 3/2 . 
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